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The edges of graphene and graphene like systems can host localized states with evanescent wave 
function with properties radically different from those of the Dirac electrons in bulk. This happens 
in a variety of situations, that are reviewed here. First, zigzag edges host a set of localized non 
dispersive state at the Dirac energy. At half hlling, it is expected that these states are prone to 
ferromagnetic instability, causing a very interesting type of edge ferromagnetism. Second, graphene 
under the influence of external perturbations can host a variety of topological insulating phases, 
including the conventional Quantum Hall effect, the Quantum Anomalous Hall (QAH) and the 
Quantum Spin Hall phase, in all of which phases conduction can only take place through topologically 
protected edge states. Here we provide an unihed vision of the properties of all these edge states, 
examined under the light of the same one orbital tight-binding model. We consider the combined 
action of interactions, spin orbit coupling and magnetic field, which produces a wealth of different 
physical phenomena. We briefly address what has been actually observed experimentally. 
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I. INTRODUCTION 

Graphene has been the most studied material of the 
last decade. Its extraordinary electronic and mechan¬ 
ical properties came as a great surprise: the existence 
of stable two dimensional crystals was customarily dis¬ 
missed, and surfaces had been identified as the source of 
reduction of electronic mobility, due to defects and ad¬ 
sorbate trapping. The age of graphene was initiated by 
the observation of the field effect transistor^, and more 
strikingly the quantum Hall effecfP^, a phenomena that 
had only been observed in high mobility semiconductor 
heterostructure^. 

Graphene is a two dimensional lattice of carbon atoms 
that form a honeycomb lattice, that can also be described 
as a triangular lattice with a two atom basis, displayed 
with different colors in Fig 1. This makes of the graphene 
honeycomb lattice a bipartite lattice, a fact that strongly 
influences its electronic properties. The electronic prop¬ 
erties of graphene can b e de scribed in terms of a very el¬ 
egant and simple pictur^^ by means of the Dirac equa¬ 
tion. Glose to the Fermi energy, electrons in graphene 
behave as two dimensional relativistic masless particles, 
the so called Dirac electrons. The energy bands are lin¬ 
ear, E± = ±TiVF\k\, so that the three dimensional plot of 
these two dimensional bands produces the so called Dirac 
cones. The Brillouin zone associated to the honeycomb 
lattice is also hexagonal, and has a copy of these Dirac 
bands, located at the corners of the hexagon. Only two 
of these so called valleys are actually non-equivalent. As 
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a result, electrons in graphene have an additional isospin, 
the valley. 

All these properties are also expected for a wider 
class of material systems, the graphene-like materials, 
that can also be described in terms of electrons mov¬ 
ing in a honeycomb lattice with just one orbital per 
site. An incomplete list of graphene-like materials in¬ 
cludes Silicen^ Germanen^, Stanen^, Metallic or¬ 
ganic frameworlPSl^ hydrogenated Bi(lll)(i3, and artifi¬ 
cial graphene lattices!^ 

The purpose of this paper is to review what is known 
about the fate of the Dirac electrons at the edges, the 
boundaries of these otherwise endless two dimensional 
crystals. In some instances Dirac electrons simply scatter 
at the edges but, very often, graphene hosts edge states, 
i.e., states whose wave function are evanescent in the 
direction perpendicular to the edge, and itinerant in the 
parallel direction. Their energies are at, or close to, the 
Dirac point, and very often their wave functions have 
peculiar properties, such as sublattice polarization, spin 
polarization or net spin current, just to mention a few. 

Edge states are particularly important when graphene 
is driven into what nowadays are known as topological 
insulator phaseJ^^. Historically, the first example of this 
phase is associated to the Quantum Hall Effect (QHEjP, 
observed in high mobility two dimensional electron gases 
in semiconductor heterostructures. In these systems, ap¬ 
plication of a sufficiently large magnetic field produces 
a discrete spectrum of Landau levels (LL) in the bulk 
states. This leads to an insulating state when the Eermi 
energy lies in between the LL. Importantly, in that situ¬ 
ation the edges host chiral (o r unid irectional) states that 
are ideal quantum conductor j^^b6 ( and are responsible of 
the perfect quantization of the Hall conductance. 

Soon after the discovery of the QHE, it was shown by 
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FIG. 1: (a) Honeycomb lattice showing two types of edge, 
zigzag and armchair. The two triangular sublattices, A and B, 
are displayed with fake color, red and blue. The vectors of the 
Bravais lattice are also shown, (b) Brillouin zone associated 
to the honeycomb lattice, including the plot of the two energy 
bands forming Dirac cones in the neighborhood of K and K' 
points (see text) 

Thouless and coworkers (TKNlSfi^ that the Hall conduc¬ 
tance could be expressed, using the conventional linear 
response theory, in terms of a topological invarianl! ^^b8( 
the so called Chern number C, associated to the Berry 
curvature of wave functions of the bulk states. 

The prediction of other insulating phases with quan¬ 
tized edge transport due to topological order without a 
net magnetic field is one of the greatest successes of mod¬ 
ern condensed matter theory. This includes the quantum 
spin Hall (QSH)p2l2ni^ quantum Anomalous Hall 

(QAH) phases, proposed in a seminal paper by Haldan^^ 
where he showed how spinless fermions moving in a hon¬ 
eycomb lattice exposed to a periodic magnetic field with 
no net flux would display quantized Hall conductance, 
with topologically protected edge states. 

The QSH phase was proposed by Kane and Mel d^^ * ^*^ . 
They found that intrinsic spin-orbit coupling would open 
a gap in graphene with non-trivial topological order that 
would come accompanied by spin filterecP^ edge states 
robust with respect to time reversal perturbations. In¬ 
terestingly, the description of electrons with spin-orbit 
coupling in graphene was mathematically identical to two 
independent copies of the Haldane model, one per spin. 
They also introduced a Z 2 topological classificatiorP^ of 
time-reversal invariant two dimensional systems, analo¬ 
gous to the TKNN classification of quantum Hall states. 

Subsequent computational worlP^HSI showed that the 
magnitude of the intrinsic spin-orbit coupling in graphene 
was so small that would render the observation of the 
QSH phase almost impossible. However, there are 
graphene-like materials, such as Silicene and other group 
IV honeycomb crystals, for which the Kane Mele model 
applies^Sl and for which these predictions are relevant. 
More importantly, there is quite strong experimental ev¬ 
idence that the QSH phase has been observed both in 
HgTe quantum wella^ and inverted InAs/GaSb quan¬ 
tum well j^, both theoretically predicted to be QSH 
insulatord^^® 

The role of Coulomb interactions can also affect dra¬ 
matically the properties of some of these edge states, 
in particular, whenever edge states produce a large 


density of states at the Fermi energy, that makes 
them prone to Stoner instabilities. This is the case 
of zigzag edge states for which ferromagnetic order 
is expectecpSHSU The interplay between this mag¬ 
netism, spin-orbit interactiond^SMSl Quantum 

Hall phase^^ is a very fascinating area of research that 
we also review here. 

Appart from the previous examples, interfacial effects 
can also create topologically protected states. Some 
examples are d riven by domain boundaries betwe en 
gapped graphen^^^EIl, local electric edge field^^^HZI or 
interfaces between antiferromagnetic graphene and a 
superconductor.!^ 

The rest of this review is organized as follows. In sec¬ 
tion H we review the tight-binding model that describes 
both the 2D and edge states in graphene, including the 
spin-orbit coupling and coupling to the magnetic field, re¬ 
sponsible of the Quantum Spin Hall and Quantum Hall 
phases. In section HI we review the properties of the 
zigzag edge states, including their connection with the 
bipartite character of the honeycomb lattice as well as 
the ferromagnetic order associated to Coulomb interac¬ 
tions. In section IV,V and VI we review the quantum 
Hall, quantum anomalous Hall and quantum spin Hall 
edge states, respectively. In sections VH and VHI we 
review the effect of Coulomb interactions on the spin- 
filtered edge states in graphene. In section IX we briefly 
review the experimental situation and in section X we 
wrap up with some general conclusions. 

II. TIGHT BINDING MODEL FOR 

GRAPHENE AND GRAPHENE-LIKE 
MATERIALS 

A material is said to be graphene-like if its quantum 
states can be described in terms of a tight-binding model 
that describes electrons hopping in a honeycomb lattice 
with a single state per site. In the case of graphene, Sil¬ 
icene, etc, the site would be a group IV atom, and the 
state would be pz orbital. Within this model, electrons 
can hop to their first neighbor atoms, with a hopping am¬ 
plitude t, that takes a value of t ~ 2.7e\® in the case of 
graphene. The tight-binding approach can also include 
the effect of magnetic fields by means of the so called 
Peierls substitution. Basically, the effect of the magnetic 
field consists on multiplying by a phase the hopping in¬ 
tegrals —>■ where 



and A is the vector potential applied to the system. 

A. Bloch and Dirac Hamiltonians 

The honeycomb lattice of bulk graphene can be treated 
as two interpenetrating triangular lattices, that we label 
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as A and B and assign them the red and blue color in 
Fig. 1. Thus, the honeycomb lattice is a triangular lat¬ 
tice with two atoms per unit cell that naturally leads, 
within the simple TB model, to a Bloch Hamiltonian of 
dimension two: 


Boik) 


t if(k) ] 

tnk) -f; 


( 2 ) 


where t is the first neighbor hopping, A is the so called 
mass term that is present whenever there is a sub¬ 
lattice symmetry breaking perturbation and it is as¬ 
sumed to vanish in the case of freestanding graphene, 
f{k) = 1-1- -I- is the form factor associated 

to first neighbor hopping in the honeycomb lattice, and 
oi = I (-x/S, l), a 2 = I (■\/3, —l) where a is the unit 
cell spacing, which coincides with the second neighbor 
distance and it satisfies the relation a = a/Succ with 
the first neighbor distance. The resulting energy bands. 




|t/(fcp are shown in Fig. 


1 D for the 


A = 0 case relevant for graphene , and feature the so 
called Dirac cones at the corners of the hexagonal Bril- 
louin Zone. Valence and conduction band meet at the 
so called Dirac point, which coincides with the Fermi 
energy at half filling. At this point we introduce the con¬ 
cept of sublattice as a pseudo spin degree of freedom. For 
that matter, we can write down the Bloch Hamiltonian 
in terms of the Pauli matrices (Ta,, cTy, Uz: 


Boik) 


t Re{f{k))a^+Im{f{k))ay 



h ■ (t(3) 


This notation makes it apparent that the Bloch states of 
graphene, that we can describe as two component spinors 

^. Actually, for A = 0, those states are analogous 

to the eigenstates of a pseudo spin under the influence of 
an in-plane magnetic field. 

As a result, the projection of their wave functions over 
the two sublattices have the same weight \(j)A\ = \4>b\- In 
other words, the bulk states of graphene are sub-lattice 
unpolarized. 

In the neighborhood of the K and K' points, denoted 
by the label = ±1, it is convenient to Taylor ex¬ 
pand the Bloch Hamiltonian to obtain the so called Dirac 
Hamiltonian: 

Bo{^ = hvF {qxCTx+ T:,qyay) + (4) 



where q = k — and vp = 3tacc/2?i, Many of the 
low energy properties of graphene can be understood in 
terms of the previous continuous model. 


B. Spin dependent terms 


Hamiltonian of electrons in graphene can have up to three 
different spin-dependent terms. At finite held, the Zee- 
man term 

Bz = ^ 91 ^bB ■ a (5) 

where g — 2. 

The other term is spin orbit coupling (SOC). The orig¬ 
inal atomic spin-orbit term, XS ■ L, has a vanishing value 
on the Pz- However, higher order processes involving or¬ 
bitals from the Px and Py manifold, or even from the d 
manifold,!^ will add an effective SOC term to the Hamil¬ 
tonian. 

Whereas a constructive procedure has not been derived 
in general, the following Hamiltonian postulated by Kane 
and Mele, reproduces all the important propertie s that 
the actual SOC brings in more realistic calculation d^^ l ^^ l 

HkM = ^ itKMO-’^a,0cl^a'^l3cr ( 6 ) 

{{a,P)),cr 

double angle brackets denote second neighbors summa¬ 
tion, a = ±1 are the spin projections (along the axis 
perpendicular to the crystal plane) and Vap = +(~)1 for 
clockwise (anticlockwise) second neighbor hopping. 

When added to the first-neighbor hopping Hamilto¬ 
nian, the Kane-Mele term opens a band-gap Asoc = 
Gy/StKM at the Dirac points. Thus, SOC would turn 
graphene into a gapped material and, as it was also found 
out by Kane and Mele, of a very special nature. In con¬ 
trast with most of the band insulators, graphene was pre¬ 
dicted to be a quantum Spin Hall insulator, i.e., topolog¬ 
ically different from vacuum. A dramatic consequence 
of the topological non-triviality of the QSH phase is the 
existence of chiral edge states in graphene. 

At first sight, it is somewhat surprising that in the 
Kane-Mele spin-orbit coupling Hamiltonian, Sz is a good 
quantum number, [Bkm,Sz] = 0. Mirror symmetry is 
the ultimate cause of this conservation law, which can 
also be understood as follows. Spin-flip terms S^L~ + 
S~L~^ connect the tt orbitals with the Px,y orbitals. Thus, 
the effective Hamiltonian has to include these processes 
in couples, so that the tt electron with spin f couples 
to a state p^ + ipy with spin j, and then return to the 
original state, preserving the spin thereby. In addition, 
density function theory calculations showfl^ that even 
though the low energy properties are dominated by the 
Pz orbitals, the SOC contribution comes also from the 
d manifold. Actually, this contribution turns out to be 
even bigger than the one from the p manifold. 

Whenever mirror symmetry is broken, due to applica¬ 
tion of an external off-plane electric field, or due to in¬ 
teraction with the substrate, another spi n orb it coupling 
arises known as the Rashba HamiltoniarP^I^ 

BR = itR E- (fij X c\^cjs' (7) 

i,j,s,s' 


The spin plays a crucial role in most of the edge state 
physics in graphene and graphene-like systems. The 


where r^- is unit vector along the bond between the car¬ 
bon sites i and j, t? are the spin Pauli matrices and E 
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is a vector related to inversion symmetry breaking of the 
graphene lattice, such as an off-plane electric fielcP^. The 
Rashba spin orbit coupling does not commute with S^. 

C. Coulomb interaction 

In general, electron-electron interactions play a sec¬ 
ondary role in the electronic properties of graphene and 
graphene-like system^^, whose defining properties are 
captured by the tight-binding model presented above. 
However, in those instances where the single-particle 
spectrum has a large degeneracy, such as the case of 
zigzag edge states as well as the bulk LL, interactions 
can have a strong effect. In order to model electron- 
electron interactions, two approximations are often em¬ 
ployed. First, only the intra-atomic part of the interac¬ 
tion is considered. This is the so called Hubbard approx¬ 
imation: 


Graphene ribbons are interesting by their own sake, and 
there has been enormous progre ss in the fabrication of 
ribbons with smooth edges^^"^, s o that inter-edge cou- 
pling is also an interesting 

A second strategy to calculate edge states is to calcu¬ 
late the Green’s function of a semi-infinite two dimen¬ 
sional crystal using the recursion method. Using the 
translation invariance along the direction parallel to the 
edge, it is possible to write the Hamiltonian of the semi¬ 
infinite crystal as a one dimensional semi-infinite crystal 
whose effective Hamiltonian depends on the transverse 
wave vector k. By so doing, the conventional Green’s 
function technique^^ normally used to deal with ID 
problems can be used. In particular, the Green’s function 
of the unit cell in the semi-infinite crystal reads: 

G{k,E) = {E - ho{k)+ie) (12) 

where S(fc) is the self-energy induced by the coupling to 
the rest of the crystal: 


T~l-u = U ^ ( 8 ) 

i 

where the sum runs over the sites i of the lattice. The re¬ 
sulting Hubbard model can not be solved exactly, except 
in a monostrand one dimensional chain. Thus, very of¬ 
ten the model is treated at the mean field approximation, 
where the exact Hamiltonian is replaced by an effective 
Hamiltonian 


'^MF ^Hartree T ^Fock (9) 

where 

^Hartree — U (f9) 


■Hpock = -U -f (II) 

so that electrons interact with an external field 
that is self-consistently calculated. In most of 
paper j^ ^° I ^^F5 | 37 | 38 | 40 | 5i | 52 ] additional approximation 
has been used, that assumes a collinear magnetization 
so that the Fock term vanishes. However, instances in 
which this is not the case are very interesting and also 
reviewed her^. 


D. Calculation method for edge states 

Two methods are normally used to compute edge 
states starting from a tight-binding Hamiltonian. In 
most instances we compute the energy bands of a so 
called graphene ribbon, a one dimensional crystal with 
two edges^. When the width of the ribbon is sufficiently 
large, compared to the penetration length of the edge 
states, the interactions between the edges are negligi¬ 
ble, allowing to study the properties of the edge states. 


Eik, E) = t{k)G{k, E)t^k) 


(13) 


where h^ik) and t{k) are the intracell and intercell ma¬ 
trices of the Block Hamiltonian. Equations (12) and (13) 


define a set of non-linear coupled matrix equations that is 
solved by numerical iteration. Once the surface Green’s 
function is derived, the density of states can be obtained 
through p(fc, E) = — ^Im {G{k, E)). Gontour plots in the 
k, E plane can reveal the existence of in-gap edge states, 
as shown below. 

Most of the results presented in this manuscript re¬ 
produce existing results of the literature. However, we 
have re-computed most of them using a home-made code. 
Quantum Honeycomp, a package that computes the 
electronic properties of honeycomb tight-binding models 
in one dimension, including the couplings to the magnetic 
field (eqs. the spin-orbit couplings (eqs. ®) as 

well as the mean field Hubbard terms in the non-collinear 
approximation (eqs. ( fiTO )■ The code is available 
onlin^^ and has a graphical user interface that permits 
a simple use by non-experts. 


III. ZIGZAG EDGE STATES 

We now review the properties of the simplest type 
of edge state in graphene, the so called zigzag edge 
states. Figure la shows the two simplest types of edges in 
graphene, the so called armchair and zigzag. Of course, 
other terminations are possible in principl^^, but we 
shall not discuss them here. As we show now, the zigzag 
terminations host E = 0 energy states, whereas the arm¬ 
chair terminations do not. The ultimate reason for these 
states is sublattice polarization. Firstly, the zigzag edges 
present a very strong sublattice imbalance. Secondly all 
the atoms of a given zigzag edge belong to the same sub¬ 
lattice. This turns out to be a very important property 
that leads to the existence of if = 0 evanescent states. 
This can be seen in several ways. 
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A. Effective mass description of Zigzag edge states 


A quick and simple argument for the existence of mid¬ 
gap E = 0 edge states in graphene can be obtained 
using the effective mass Hamiltonian from eq. (18) 
taking E = 0 and B = 0 and proposing the antsaz 

which the wave function lives only in 




one sublattice. It is apparent that <^o{x) = satisfies 
the Dirac equation. Therefore, the localization length 
can be written as A = q~^, is maximal for the Dirac 
point, g = 0, and decreases (the state becomes more lo¬ 
calized) as q increases. A more complete analysis along 
this line, including the calculation of the inter edge cou¬ 
pling in the case of graphene ribbons, can be found in the 
seminal work of Brey and Ferti^^. The validity of the 
effective mass approximation to describe states abrupt 
perturbations, such as an edge, is questionable a priori. 
However the results of this approximation compare well 
with the tight-binding results^. 


miLLi 
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FIG. 2: (a) Calculated magnetic moments Zigzag edge states 
for a graphene ribbon with A = 20 atoms in the unit cell, both 
in the AF and FM configurations, using the meant field Hub¬ 
bard model with U = t. Energy bands for the non-interacting 
(b), and the AF (c) and FM (d) solutions for a ribbon with 
A = 40 atoms in the unit cell. The color red (blue) denotes 
spin up (down), notice that in figure (a) only one of the spins 
are visible but the degeneracy is two for the whole spectrum. 


B. Tight-binding description of Zigzag edge states 

The existence of an A = 0 edge state was first in¬ 
ferred from the calculation of the energy bands of a 
zigzag graphene ribbon in the early work of Nakada and 
coworkerjA^. Direct diagonalization of the TB Hamilto¬ 
nian yields the energy bands that we show in figure &)■ 
Two types of bands are seen: confined bulk modes, that 
have a gap, and two flat bands with A ~ 0 corresponding 
to states localized at the two edges of the ribbon. Small 
departures of their energy from A = 0 are related to inter 
edge hybridization. This occurs because the edge wave 
function associated to a given edge lives in the opposite 
sublattice than the edge state coming from the opposite 
edge. An important property of the electronic states of 
zigzag ribbons is that they do not mix valley^^. Actu¬ 
ally, the confined bulk modes are cuts of the two Dirac 
cones, whereas the A ~ 0 flat bands join the otherwise 
disconnected valleys. 

The A = 0 wave function of a single edge can be 
worked out analyticalljff^. If the edge atoms belong to 
the A sublattice, they found that (fiB = 0 and (j)A{rn) = 

(2cos where m = 1, 2,.. labels the A atoms start¬ 

ing from the edge along an armchair row. This wave 
function is only normalized if fca > i.e., between the 

Dirac points and the T point, which naturally explains 
the results of figure [^. An expansion of the wave func¬ 
tion coefficient around the Dirac point gives an exponen¬ 
tial decay, similar to the one obtained using the Dirac 
Hamiltonian. 


C. Zero modes in bipartite lattices 

The zero modes of the zigzag edge can be understood 
in a broader context. A one orbital tight-binding Hamil¬ 
tonian defined on any bipartite lattice, with a number of 
atoms in one lattice larger than in the other {Na > Nb), 
has least 

Nz = \Na — Nb\ (14) 

states with zero energy whose wave function is local¬ 
ized in the A (majority) sublattic^^. Three standard 
examples of systems with sublattice imbalance would be 
graphene with a vacancj^^, semi-infinite graphene with a 
zigzag edge and a triangular graphene island with zigzag 
edge^. 

The proof of this theorem goes as follows. For a bipar¬ 
tite lattice, the Schrodinger equation can be written as 
two linear systems: 

Hab‘^b = E^a 
b 

Y Hba<^a = E<^b (15) 

a 

In order to verify the theorem, we take A = 0 and (f>b = 0 
for the Nb sites of the lattice, so that we are left with an 
undetermined linear homogeneous system of equations, 
J2a Hba^a = 0, with Na variables, the components of the 
wave function on the majority sublattice, but only Nb 
equations. There are an infinity of such solutions, which 
form a vector space, whose dimension is the difference 
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between the number of unknowns Na and the rank of 
the matrix of the system, Nb- This vector space is the 
one of the E = 0 modes enunciated in the theorem. 

In the case of a bipartite lattice with a constant stag¬ 
gered potential, such that all A (B) atoms have an on-site 
energy ^ (^■) the theorem can be also applied, but now 
the ’’zero modes’ have energy ±^, depending on which 
one is the majority sublattic^^, and they are still entirely 
localized on a single sublattice. 

A dramatic example of A = 0 edge states is provided 
by triangular graphene nano islands with zigzag edge^^, 
shown in figure A remarkable feature of this class 
of systems is that their sublattice imbalance scales with 
size, and so it does the number of A = 0 states. So, the 
triangular islands can have strict E = 0 states, even for 
small systems. The smallest of such systems has an odd 
number of A^ = 13 atoms, so that it has to have a sub¬ 
lattice imbalance, Na — Nb = 1- But the second of such 
structures has an even number of atoms, N = 22, which 
would permit to build many structures with Na = Nb 
such as pentancene , yet the triangulate has Na — Nb = 2 
and, thereby two mid-gap E = 0 states. Inspection of the 
wave function reveal that their wave functions have more 
weight on the edge^^. Another example of ill = 0 edge 
state is provided by the minimal edge that can be de¬ 
vised, the one created by a single vacancy. The theorem 
warrants the existence of a £1 = 0 state. The lack of a 
gap in graphene turns this E = 0 state into a resonance 
whose amplitude decays . In the case the quantum spin 
Hall phase, where a non-trivial gap is open, the vacancy 
creates a real mid-gap state with a normalizable wave 
functiorPS. In both cases, the vacancy localizes an elec¬ 
tron 

It must be stressed that the theorem of eq. (14) gives 
the minimal number of zero modes of a given structure. 
As discussed by Koshino and coworker^, the presence 
of symmetries in the Hamiltonian that commute with 
the chiral symmetry, the one responsible of the electron- 
hole symmetry and other features of bipartite Hamiltoni¬ 
ans, can increase the number of zero modes if the invari¬ 
ant subspaces associated to that symmetry have their 
own sublattice imbalance. As a result, structures with 
a global null sublattice imbalance still have zero energy 
statePl. 


D. Edge magnetism: mean field Hubbard model 

The E = 0 states discussed in the previous subsection 
play a very important role in graphene because, at half 
filling, the Fermi energy lies exactly at A = 0. Actually, 
in a system with Na + Nb atoms at half filling, with 
Na > Nb, the number of bound states with A < 0 is 
Nb, so that the valence band hosts 2Nb- If we write 
the number of electrons N = Na + Nb = Nz + 2 Nb, 
it is apparent that the in-gap E = 0 states are half¬ 
full. Given that the wave functions of the E = 0 states 
overlap in space, the Coulomb interactions are expected 


to favor ferromagnetic spin correlations, very much like 
in open shell atoms. In a solid state physics parlance, the 
presence of edge states creates a large density of states at 
the Fermi energy. Therefore, Coulomb interactions could 
result in a Stoner instability that produces ferromagnetic 
order. 

These hand-waving arguments can be put on a very 
firm basis. Density functional theory (DFT) cal¬ 
culations have shown that zigzag edges are indeed 
ferromagneticpi^^^. In the case of graphene ribbons, 
the magnetization of the edges affects the otherwise flat 
edge bands, which become dispersive. When the mag¬ 
netization of the edges is (anti)parallel, the so called 
(anti)ferromagnetic configuration, the ribbon is (non) 
conducting. The change of the conduction properties 
of a zigzag graphene ribbon with the relative orienta¬ 
tion of the edge magnetiz ation inspired the proposal of a 
graphene based spin valv d^^ * ^^ as well as a Silicene spin 
valv^^. 

Interestingly, the description of the magnetic order in 
graphene ribbon u sing t he mean field approximation of 
the Hubbard modePSM] gave results very similar to those 
of DFT and, in addition, provided analytic insight on 
origin of the magnetization induced dispersion as well 
as a treatable description of the evolution of magnetism 
away from half fillin^^ . The difference between the 
U — 0 bands (fig. (|^)) and the U ^ 0 bands (fig. 
([^,d)) is known as the interaction self-energy. In the 
case of the AF configuration, the self-energy is non zero 
only for the two edge bands, so that it can be calculated 
analyticalljEH. In that case the self-energy is a dimension 
two matrix and it has both intra-edge and inter-edge con¬ 
tributions. The former gives rise to the spin splitting of 
the band states close to the F point k = tt and is thereby 
identical for the AF and FM solutions. As we move away 
from the F point towards the valleys, the wave functions 
become more delocalized and the inter-edge contribution 
to the self-energy takes over and, depending on the rel¬ 
ative orientation of the edge magnetic moments a gap 
opens (AF) or a spin polarized conducting channel re¬ 
mains open at each valley (FM). 

The results of the mean field Hubbard model calcula¬ 
tions for the ferromagnetic (FM) and antiferromagnetic 
(FM) for a graphene ribbon are shown in figure tak¬ 
ing {7 = t for a ribbon with A^ = 40 atoms. For this 
value oi U/t = 1 the magnetization of the edge atoms is 
0.15 ^b, three times smaller than the magnetic moment 
of Nickel atoms in ferromagnetic Nickel. The magne¬ 
tization of non-edge atoms is much smaller and decays 
exponentially as a function to the distance to the edge. 
Within the same mean field Hubbard model, it would 
take a value of U larger than 2.2t, to produce magnetic 
order in the bulk honeycomb lattice, which would be anti¬ 
ferromagnetic (AF). More sophisticated approximations, 
such as quantum Monte Carlo, push the critical value of 
U to even higher value^^. In the case of edge atoms, 
the mean field magnetization is non-zero for any positive 
value of L®. 
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FIG. 3: (a) Non-interacting {U = 0) energy levels for the 
N = 22 triangulene, Notice the two states with E = 0. 
(b) Calculated magnetization, within the mean field Hubbard 
model with U = t. The self-consistent solution has two un¬ 
paired electrons {S^ = 2) whose magnetization goes predomi¬ 
nantly to the majority sublattice atoms (shown in red) at the 
edges. 

In the case of triangular graphene islands, both density 
functional calculations and mean field Hubbard model 
calculations yield very similar result^^. Triangular is¬ 
lands whose single particle spectrum has Nz zero en¬ 
ergy states develop edge magnetism with Sz = as 
if these systems were following the atomic Hund rul^^. 
Actually, in the case of the Hubbard model, Lieb demon¬ 
strated the following theorem: the spin S of the ground 
state of a Hubbard model defined in a bipartite lattice, 
at half filling, is given by 5" = AavAb , Thus, the mean 
field approximation is, in some sense, compatible with 
the Lieb theorem, but it is important to keep in mind 
the differences between the broken symmetry mean field 
solutions and the exact solutions. In the case of graphene 
ribbons, the AF configuration has lower energy than the 
FM solution, which very often is said to be in agreement 
with the Lieb theorem, but the exact solution with S' = 0 
would have a vanishing magnetic moment at each edge, 
whereas the mean field AF solution the net magnetiza¬ 
tion vanishes, but not the local magnetization. 


E. Beyond the mean field Hubbard model 

The mean field theory of the Hubbard model provides 
a handy first description of the magnetic properties of the 
edges of graphene and graphene-like materials. However, 
this approach has several well known shortcomings. For 
instance, the long-range tail of Coulomb interactions in 
graphene should be poorly screened. Therefore, non-local 
exchange effects are missed both by mean field Hubbard 
approximation and DFT calculations based on local den¬ 
sity approximations. A proper treatment of non-local 
exchang^^ shows that magnetic features are further sta¬ 
bilized, both for intra-edge and inter edge exchange. 

Both local and non-local mean field approximations 
ignore spin dynamics and the fact that long-range or¬ 
der is not possible in one dimension. The spin waves of 
the magnetically ordered zigzag ribbon were calculated 
within the RPA approximation for the Hubbard model by 
Wakabayashi and coworker^^ who found gap-less Gold- 


stone modes, as expected. In consequence, long-range 
order should be suppressed. Therefore, at finite tempera¬ 
ture intra-edge spins correlations are expected to survive 
up to a temperature dependent spin correlation length 
that was computed by Yazyev and Katsnelsorf^. At 
room temperature the magnetic correlation length was 
estimated to be Inm. 

Within the mean field approximation the Fermi liquid 
picture of quasiparticles remains intact. However, in one 
dimension, interactions are expected to fully renormalize 
the one-particle bands, and strongly correlated phenom¬ 
ena such as spin-charge separation are expected to occur. 
This scenario has been considered for the case of quantum 
Hall edge state^^, as well as the quantum Spin Hall edge 
stateP^ and more recently for the case of ferromagnetic 
zigzag edges as welP^. 


IV. QUANTUM HALL PHASE 

In this section we discuss a totally different type of 
edge states. They arise in graphene when a magnetic 
field is applied perpendicular to its surface, they are ex¬ 
tended much more than a few unit cells, and they can 
live in all types of edges, not only at zigzag boundaries. 
These edge states are associated to the so called Quan¬ 
tum Hall state, which shows a bulk insulating behavior 
together with perfectly conducting edge transport. In a 
classical naive picture, bulk insulating behavior is asso¬ 
ciated to electrons performing closed orbits, whereas the 
edge electrons are able to move forward in the so called 
skipping orbits by bouncing on the interface. 

Nevertheless, the quantum mechanical treatment gives 
a much richer understanding of QHE. Bulk localization 
is associated to the emergence of the so called Landau 
levels (LL), with a discrete spectrum that replaces the 
energy bands and localized wave functions. The perfectly 
quantized edge conductance can be understood by taking 
into account the non-trivial topology of the bulk states 
as well as the bulk to edge correspondence. 

The role of topology can be seen as follows. Both vac¬ 
uum, or of that matter any trivial insulator, and the 
quantum Hall state are insulators. However, there is 
no way to change the parameters in their Hamiltoni¬ 
ans to turn one into the other without closing the band 
gap. Therefore, the space of Hamiltonians generated by 
changing parameters in the trivial and quantum Hall in¬ 
sulators, without closing the gap, gives rise to (at least) 
two disconnected classes of insulating Hamiltonians that 
can not be deformed one into each other without pass¬ 
ing through a gapless state. As a consequence, at the 
physical interface between two materials described with 
Hamiltonians that belong to these two different classes 
there must be a conducting state. The general argument 
does not only apply to the quantum Hall state, but also 
to a large amount of different systems known as topologi¬ 
cal insulators. In the next sections of this review, we shall 
discuss also some of them in the context of graphene. 








(17) 


A. Electronic states 


valley!^ 


The description of the electronic states of graphene 
in a constant magnetic field starts with the choice of a 
vector potential A = B(—y,0,0) to be inserted in eq. 
Q, so that the original periodicity of the crystal is only 
preserved along the x direction. As a result, the TB 
description of graphene under a magnetic field is very 
often restricted to one dimensional stripes which permits 
to study the bulk and edge states on equal footing. In 
a tight binding model this minimal coupling is realized 
via Peierls substitution. The magnetic field introduces a 
new length scale in the problem: 



which is Ib = 25.7nm for B = IT. For a sufficiently high 
magnetic field the magnetic length is much smaller than 
the width of the system W. In that situation, the bulk 
quantum states become localized, giving rise to LL, but 
the edge states are dispersive, as shown in figure ob¬ 
tained by numerical diagonalization of the TB Hamilto¬ 
nian for a ID stripe with zigzag edges (left) and armchair 
edges (right). 

The spread of the edge wave functions depends on their 
energy with respect to their reference Landau leveP^, 
scales with Ib and is thereby much larger than the spread 
of the B — 0 zigzag edge states discussed in the previous 
section. Importantly, in the case of zigzag edges, both 
the magnetic edge states and the E = 0 edge states co¬ 
exist. Actually, inspection of the wave functions show 
that, at a given valley, the two bands with E = 0 bands 
shown in fig.([^) correspond both to the n = 0 Landau 
level in that valley and the edge state. 


0.4 

0.2 

E/to.o 

- 0.2 
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FIG. 4: LL for zigzag (a) and armchair (b) ribbons. Band 
structure for a graphene ribbon with orbital magnetic field. 
The color code means red and blue for the edges and green 
denotes the bulk states. 

Further insight of these results can be obtained within 
the so called effective mass description, that can be for¬ 
mally derived from the kp theory of the energy bands. 
This technique had tremendous importance in the de¬ 
scription of semiconductor^^. In practical matters, it 
amounts to replace q in equation by the momen¬ 
tum operator. By so doing, we obtain an effective mass 
Hamiltonian isomorphic to the Dirac Hamiltonian at each 



0.2 0.4 0.6 0.8 

fca/( 27 r) 



0,2 0.4 0.6 0.8 

ka/{ 27 T) 


Hr = Vp + TUyay) -I- Acr^, 

where H = p — eA is the canonical momentum opera¬ 
tor, vp = 3tacc/‘^T^i ^ are the Pauli matrices describing 
the graphene sublattice degree of freedom and r = ±1 
describes the valley index. 

In the case of constant magnetic held perpendicular 
to the plane this Hamiltonian can be solved exactly. Us¬ 
ing the transnational invariance along the x direction, we 
can assume its eigenfunctions are products e^’^^^(j)n{kx,y) 
which permits replacing the operator by the quantum 
number fikx in Eq. (Hii- The Schrodinger equation be¬ 
comes then a problem of two Hrst order differential equa¬ 
tions in one dimension: 

Hr = vpifikx - eBy)ax + rvppyay) + (18) 

The solution of the problem is facilitated by introduc¬ 
ing the dimensionless operators: Q{kx) = — kxlp^ 

and P = ^-§Py, and the ladder operators 

a{kx) =-^{Qikx)+ iP), (19) 

which satisfy bosonic commutation relations 

[a{kx),a{kx)^] _ 1. Using these operators, the 

Hamiltonian ( |l^ for valley r = — 1 can be written as: 

H = ^ (cr+a -b cr"a1') (20) 

where = (Ja, ± iCy where we have dehned = 

For T = -bl we have to exchange a and ab 
Tfiis model is mathematically equivalent to the Jaynes- 
Cummings model in quantum optics describing a two 
level system coupled to a bosonic mode. The result¬ 
ing eigen function s and eigenvalues have the following 
propertie d"^^ l ^^ l ^^ [ First, there is a spectrum of states 
with wave functions with non-zero (j)A and with en¬ 
ergies: 


EN = ±\lA^ + ^{huJofN 


( 21 ) 


with N a strictly positive integer. This spectrum comes 
with a twofold valley degeneracy, in addition to the 
twofold spin degeneracy. Notice, that for A = 0 the Lan- 
dau level spectrum for Dirac particles scales with ^/NB, 
in contrast to the linear scaling of Schrodinger particles, 
that can be retrieved in the limit of very large A. 

In addition to the states described with eq. (21), there 
is one zero mode per valley, with energy Eq = and 
with wave function fully sublattice polarized: 


[>T = + 1 — 


4>o 

0 


»r=-l — 


0 

4>o 


( 22 ) 
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For A = 0, these zero mode would be degenerate at E = 
0. Importantly, the bulk Dirac-Landau levels obtained 
analytically from the effective mass approach are in very 
good agreement with the tight-binding calculation. For 
the edge states it is also possible to work out the kp 
theory^, but this goes beyond the scope of this review. 


B. Topological origin of the edge states 

As we anticipated at the beginning of this section, in 
the context of the non trivial topology of the band struc¬ 
ture, the edge currents are just protected interface states 
in a boundary between two insulators with topologically 
inequivalent ground states. Moreover, the perfect con¬ 
ductance of this state, is a result of the relation between 
the quantum Hall conductance and the so called Chern 
number 

C^xy = j^-C (23) 

where C is the flux over Brillouin zone 

C = ^ [ ndPk (24) 

27r J 

of the Berry curvature H 

H (25) 

which acts as an anomalous term in the velocity. 

The Chern number of normal insulators is C = 0. In 
contrast, the Chern number of the quantum Hall state is 
an integer different from zero, which results in the quan¬ 
tization of the Hall conductance. 

In the case of interfaces between two different quantum 
Hall states, the previous interpretation gives a simple way 
to predict the number of edge states. For an interface be¬ 
tween a first system with Chern number C\ and a second 
system with Chern number C 2 , the number of protected 
edge states is simply |Ci — C 2 I. 

In the quantum Hall regime, the lack of edge states 
when a subl attice imbalance mass gap opens at half- 
fillinJ^HZlIZI!^ accounts for C = 0 and can be rational¬ 
ized in terms of adiabatic evolution of the massive Dirac 
Hamiltonian from B = 0. Taking the half-filled case as 
a reference, the total Chern number away from half fill¬ 
ing can be calculated by adding one unit for each Dirac 
Landau level crossed by the Fermi energy as we move it 
up, including the valley and spin degeneracy. Conversely, 
if we consider hole doping, the Chern number becomes 
negative as we move the Fermi energy down in energy. 
For example, starting from half filling filling the electron 
n = 0 Landau level gives C = 1 per spin channel, as can 
be observed in Fig. 4. Analogously, un-filling the hole 
n = 0 level gives a total C = — 1, per spin channel. This 
will lead to remarkable phenomena when dealing with 
spin polarized systems. As we discuss below, different 
fillings for the up and down channels, can give rise to a 
very special type of quantum spin Hall effecfP^. 


V. QUANTUM ANOMALOUS HALL PHASE 

Another kind of topologically non-trivial quantized 
phase is the so called Quantum Anomalous Hall (QAH) 
phase. This topologically insulating phase has properties 
analogous to those of the quantum Hall phase, including 
chiral edge states for which backscattering is impossible, 
and breaking of time reversal symmetry. However, the 
QAH phase is driven by internal degrees of freedom and 
interactions, such as spin-orbit coupling and exchange, 
rather than by the application of an external magnetic 
field. The fabrication of materials showing this phase 
has enormous potential since they are expected to show 
a perfect edge conductance as the quantum Hall state, 
but without the drawback of having to apply large mag¬ 
netic fields. 

Two important graphene related models present the 
QAH phase. First, the Haldane model for spinless 
fermions moving in a honeycomb lattic^^. Second, 
graphene perturbed both with a Rashba spin-orbit in¬ 
teraction and subject to an ex chang e field that couples 
to the off-plane spin componenlPSEII, 
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FIG. 5: Quantum Anomalous Hall phase in two toy models. 
Left panels: Haldane modeP^. Right panels: model proposed 
by Qiao and coworkerPanels (a,b): Bulk density of states 
showing a gapped spectrum resolved in the Panels (c,d): 
band structure of an armchair ribbon system for the same two 
models, showing that in addition to the gapped bulk states, 
gap-less edge states show up. Panels (e,f): density of states 
in the armchair edge of a semi-infinite plane, showing a coex¬ 
istence of bulk and one-way edge modes, calculated using the 
Green’s function approach. 
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A. Haldane model for QAH 

The toy model proposed by Haldan^^ consists of a 
honeycomb lattice in which there are local magnetic fields 
of different sign within the hexagon, but whose flux over 
it is zero. This gives rise to an imaginary second neighbor 
hopping, which accounts for the local magnetic fields: 

'Hqahi = Ho + t2 ^ ivijclcj (26) 

((id)) 

where = (^di x d 2 ^ ■ z is the chirality of the path of 

the second neighbor hopping, where di ^2 are defined as 
follows: for a given pair of second neighbor atoms i and j, 
with a common first neighbor k, di ^2 are the unit vectors 
along the bonds ik, kj, respectively. 

In this model, the bulk is characterized by a gapped 
spectrum as shown in (Figj^). A kp expansion of the 
t 2 term around the Dirac point would yield the following 
extra term in the Dirac Hamiltonian: 

SH = ^V3t2Tza2 (27) 

where the mass can be written as a valley draendent 
object, Aqahi = 3\/3t 2, in contrast with eq. W. This 
has a very important consmuence. The Chern number 
associated to Hamiltonian Q read^^ 

C(A,t^) = (28) 

Therefore, in the Haldane model both valleys contribute 
with a Chern number with the same sign that sum C = 1, 
whereas in graphene with a trivial gap, opened with a 
term that breaks sublattice symmetry, the two valleys 
contribute with a Chern number with opposite sign, that 
give C = 0. 

From the C = 1 for the Haldane model we expect the 
existence of an edge state. This is seen in the spectrum 
of a ribbon with armchair edges as in-gap chiral state 
in (Fig[^), together with the confined modes (in green). 
The chiral character of the edge states prevents backscat- 
tering. Therefore, the Haldane model describes a phe¬ 
nomenology identical to the quantum Hall phase, but 
without the presence of flat Landau levels in the bulk. In 
addition, a generalization of the Haldane model played a 
crucial role in the p ropos al in the proposal of the quan¬ 
tum Spin Hall phas^^^^, as we discuss below. 


B. Exchange plus Rashba model for QAH 

A more realistic model describing another realization 
of the quantum anomalous H all sta te in graphene was 
proposed by Qiao and coworker^^SM] jn their model they 
describe electrons in graphene including both the Rashba 


spin-orbit term (eq. 0 ) and exchange fields that splits 
the spin along the off-plane direction: 

Hqah 2 = Ho + Hr + X''^cl^SzCis (29) 


A qualitative understanding of how these two terms open 
a gap in the bulk structure, can be gained from the fol¬ 
lowing argume nt. T he off-plane exchange field will split 
the Dirac coned^SlSil^ go that the f and | cones are degen¬ 
erate on a circle in momentum space that, at half filling, 
corresponds to the Fermi circle. The Rashba term cou¬ 
ples t and I, opening a band gap in the Fermi circle. 


The k^ resolved bulk DOS shows again a gapped spec¬ 
trum (Figj^). However for a finite ribbon, several bands 
crossing the Fermi energy appear (Figj^). From those 
four bands that appear in the ribbon, two of them are lo¬ 
cated in one edge of the ribbon and two in the other. This 
becomes clear calculating the DOS on the edges, show¬ 
ing also that the two bands present are co-propagating, 
thus yielding a protection against back-scattering as in 
the Haldane model. The Chern number of this is system 
is C = 2,1^ in agreement with the two states per edge 
observed in Fig. il. 


There are several proposals to modify graphene so that 
it is described by Hamiltonian (m. (29)), including the 
use of ferromagnetic substrates,®! magnetic ad-atoms,®! 
and even antiferromagnetic substratea^Sl. Finally, a sim¬ 
ilar QAH can be obtained in graphen^^ by replacing the 
Rashba-like term by a topological spin texture known as 
skyrmion. 


VI. QUANTUM SPIN HALL PHASE 


The predictiorP^lsnilTl and discover^®! of the Quantum 
Spin Hall phase boosted the research field of topological 
insulators. Very much like the QH and QAH phases, 
the QSH phase has an insulating bulk and conduct¬ 
ing chiral edge states which, in contrast with the QH 
and QAH case, have well defined spin chirality: a given 
edge hosts two counterpropataging states with opposite 
spin projec tion, which has potential for spin electronics 
application®®!. Here we discuss two of the pioneer pro¬ 
posals for the QSH, both of them base d on g raphene. In 
the first one, Kane and Mele proposec®®! that intrin¬ 
sic spin-orbit coupling would turn graphene into QSH. 
Unlike the QH and QAH, this system would be time re¬ 
versal invariant. The second proposal, by Abanin and 
coworker®! is rather special: ferromagnetic order of the 
QH phase in graphene at half filling, for which C = 0, 
would results in a QSH phase, with = —C^ = -fl, but 
with time reversal symmetry clearly broken. We now re¬ 
view the most remarkable features of the edge states of 
these two fascinating QSH states. 
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A. QSH driven by spin-orbit 


The influence of atomic spin-orbit coupling on the 
Dirac bands, formed by = 0 atomic orbitals, is sub¬ 
tle: to lowest order, the atomic spin orbit coupling XL ■ S 
within the p-manifold has no effect at all. Second order 
interband scattering gives the first non-zero contribution 
but it is not obvious how to include it in the Hamiltonian. 
Further work,l^ showed that inclusion of d — channels 
lead to a first order contribution in SOC. In their semi¬ 
nal work, based on symmetry considerations, Kane and 
Mele postulated that the effective Hamiltonian for SOC 
in the subspace of the tt orbitals would be given by a spin 
dependent second neighbor hopping, eq. § , which turns 
out to be mathematically identical to the Haldane term 
that, as discussed in the previous section, turns spin¬ 
less fermions into the QAH phase. Thus, the Kane-Mele 
Hamiltonian Hq + Hkm for graphene with SOC is made 
of two decoupled copies of the Haldane model, one per 
spin, with opposite sign of the gap-opening term. Thus, 
the QSH phase described by the Kane-Mele model can 
be thought of as two QAH, one per spin, with opposite 
magnetization. 

In the neighborhood of the Dirac points, the KM model 
would give the following Hamiltonian: 


= 'Ho{q) + -V^tKMTS^a^Sz 


(30) 


where 'Ho{q) is given by eq. Q. Thus, the Kane-Mele 
gap-opening term respects both time-reversal symmetry 
and inversion symmetry, in contrast with eq. (27). 


The fact that the Kane-Mele model consists of two 
copies of the Haldane model permits to anticipate its 
most salient electronic properties: bulk states present a 
gap while the edges host one chiral conducting states 
per spin channel. Since the sign of the gap is given by 
the projection of the spin Sz, spin f and j, edge states 
propagate in opposite directions, i.e., they have a well 
defined spin chirality. As a result, the edge states carry 
an equilibrium spin current, in contrast with the QH and 
QAH phases for which edge states carry a net charge 
current (see figs. |^a,b,d)). 

Trivially, the Chern number of each spin channel is 
Cs^ = sgn{Sz), so that the spin Chern number C^—Ci = 2 
would be finite but the the total Chern number would 
vanish. So, it would seem that the addition of spin-flip 
perturbations, such as the Rashba Hamiltonian, would 
result in intra-edge backscattering. However, this is not 
the case and it turns out that the Spin Hall phase is topo¬ 
logically different from the normal phase, even when Sz is 
not conserved. To show this, Kane and Mele introduced 
another topological invariantP^ in order to classify these 
systems, the Z 2 invariant. The Z 2 group consists only 
of two elements, say, {0,1} and according to Kane and 
Mele’s classification, all non-trivial insulators would have 
a Z 2 index equal to 1 while trivial insulator would have 
Z 2 index equal to 0. When Z 2 = I the edge states are 
counterpropataging Kramers pairs, so that only pertur¬ 
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FIG. 6: Top panels: Electronic band structure of a graphene 
zig-zag ribbon with a small staggered potential and an ex¬ 
ternal magnetic field (off-plane) to resolve the degeneracies. 
Panel (a) shows the density of states at the edge of a semi¬ 
infinite plane with zig-zag termination. Panel (b) corresponds 
to the band structure of a zig-zag ribbon, the color represents 
spin polarization. Panel (c): the band structure for a ribbon 
in the presence of an off-plane Zeeman magnetic field and 
an orbital magnetic flux leading to the appearance of LL. In 
panel (d) a scheme of the QSH effect, where two spin polarized 
conducting states are present in each edge. 


bations that break time reversal symmetry can produce 
backscattering. Importantly, Kane and Mele showed that 
the addition of a small Rashba term to their Hamiltonian 
would keep Z 2 = 1, and only when the Rashba term is 
large enough as to close the bulk gap would the system 
turn into a. Z 2 = Q trivial phase without topologically 
protected edge states. 

Whereas the observation of the spin-orbit coupling 
driven QSH phase in graphene is very difficult due to 
the small size of the spin-orbit gap (smaller than O.I 
meV^ , the observation of the Quantum Spin Hall phase 
HgTe/CdTe quantum well^^ as well as in InAs/GaSb 
inverted quantum well^^ has confirmed the existence of 
this fascinating class of materials. 


B. Quantum spin Hall effect without spin-orbit 
coupling 


Even if spin-orbit coupling can not drive graphene into 
the QSH phase, there is a different mechanis m that c an 
do the trick and, contrary to other proposali^^ISnilZl^ it 
does not require spin-orbit couplinj^. This mechanism 
relies instead on the very special Landau level structure 
that emerges upon application of a magnetic field. 

Lets take graphene under the influence of an off-plane 
magnetic field. At half filling, the Chern number of all 
the occupied bands is zero and two of the four n = 0 LL, 
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valence and conduction with spin degeneracy. Due to the 
four n = 0 are degenerate, there are several alternative 
filling patterns. At half filling, all possible combinations 
give C = 0, but it is possible to have Cs = 2 if we fill 
the two LL with the same spin, and leave empty the 
other two. Interestingly, this spi n-polarized filling can 
be driven by Coulomb interactiona^SESElMl and also by 
intrinsic spin-orbit couplin^^. 

Thus, at half-filling, the spin-polarized quantum Hall 
phase has C = 0 and Cs = ±2, so that it has the same 
topological properties than the QSH phase described by 
the Kane-Mele model. The energy bands corresponding 
to the addition of a large Zeeman splitting are shown 
in figure [^c). They feature a gapped bulk together 
with counter-propagating spin-filtered in-gap edge states. 
Similar results can be obtainecP^ using a self-consistent 
calculation of the mean field Hubbard model, including 
an in-plane Zeeman field that favors ferromagnetic order. 

Unlike the QSH phase driven by spin-orbit cou¬ 
pling, spin-flip perturbations give rise to spin-flip edge 
backscattering and, because time reversal symmetry is 
broken, it is not possible to accommodate this QSH sys¬ 
tem into the Z 2 classification of Kane and Mele. Interest¬ 
ingly, there is very strong experimental evidence that this 
interaction driven QSH phase has been observed in the 
experiments, where the ferromagnetic order is favored by 
applying a quite large in-plane magnetic fielcP^. At half 
filling, the experiments show how two terminal conduc¬ 
tance of graphene is tuned from 0, corresponding to a 

2 

gapped quantum Hall phase to 1.8^, that most likely is 
the QSH phase predictecP^ by Abanin and coworkers. 


VII. COULOMB DRIVEN BREAKDOWN OF 
THE EDGE CONDUCTION IN THE QSH PHASE 

The gapless edge states in the quantum spin Hall states 
are protected by the spin Chern number in the case of 
magnetic field driven phase,^ and by the Z 2 invariant 
in the case of the SOC driven.l^ These topological in¬ 
variants are well defined as long as Sz is a good quan¬ 
tum number (Spin Chern) and time reversal symmetry is 
present (.^ 2 )- In the case spontaneous magnetism shows 
up, such conditions no longer apply, so the gap-less edge 
states are no longer guaranteed. 

In this section we show two examples of how Coulomb 
interactions can lead to the disappearance of spin chiral 
edge states due to the emergence of spontaneous edge 
magnetization. Following the structure of the previous 
section, we consider both the SOC driven, and the quan¬ 
tum Hall ferromagnet QSH phases. 


A. Kane-Mele-Hubbard model 

In order to study the interplay between edge mag¬ 
netism and topologically protected edge states we use 


the so called Kane-Mele-Hubbard modeP^ElIlSMI Al¬ 
though the gapless edge states are protected against time 
reversal perturbations, an important issue is whether 
electron-electron interactions are able to drive the sys¬ 
tem into a symmetry breaking state. The study of the 
bulk properties of the Kane-Mele-Hubbard model has at¬ 
tracted considerable attentiorPsHUl. Here we focus on 
the competition between SOC and magnetic order at the 
edgePS™. 

Edge magnetic order would break time reversal sym¬ 
metry and could gap out those edge states. E arly w ork 
assumed that magnetic order would be off-planIn 
this situation, the gapless edge states survive (Figj^), at 
least for small U and large SOCP^. In general, magnetic 
order in a topological insulator protected by time reversal 
symmetry should destroy the gapless edge states. This 
case is exceptional because in addition to being protected 
by time reversal symmetry, the Kane-Mele model shows 
a well defined spin Chern number, which relies only in 
conservation of Sz, preserved by the off-plane magnetic 
order. 

The situation changes r adicall y when the edge mag¬ 
netic moments lie in-plane! ^^boi | Then, Sz is no longer 
a good quantum number, and thus the protection that 
remained in the off-plane magnetic case no longer holds. 
Now, the in-plane magnetism opens up a spin mixing 
channel, driving the system into an edge insulator, as 
shown in Figj^. 

Interestingly, spin orbit coupling produces magnetic 
anisotropy, so that both configurations are not expected 
to be energetically equivalent. Actually, the ground state 
of the mean field Hubbard model is in fact the in-plane 
magnetic state.l^ However, the disappearance of the chi¬ 
ral edge states can be restored, by overcoming the mag¬ 
netic anisotropy of the system, and forcing the states to 
lie off-plane. This could be achieved by a magnetic sub¬ 
strate, which if made switchable, would allow to control 
the flow of the spin channels along the edge at will. 


B. Magnetic field driven QSH phase 

In this situation, the shifting between the LL driven 
by the magnetic field yields the effective Quantum Spin 
Hall state. 

The state arises when the in-plane field is large eno ugh 
to overcome the possible antiferromagnetic state.^^ 
However two situations of edge states are possible. On 
one hand, at weak interactions the edge becomes ferro¬ 
magnetic at the same time as the bulk. On the other 
hand, when interactions are stronger, the edge of the 
system is able to develop a canted magnetic order, even 
though the bulk remains ferromagnetic. 

Experimentally, it was observed that in order to reach 
the QSH in graphene, a large magnetic field of 20 T 
had to be applied.!^ Two different scenarios are possi¬ 
ble to understand the large field needed to overcome a 
the fully gapped phase. The first one consists on hav- 
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FIG. 7: Band structure (a) and local magnetization (b) of 
the ferromagnetic quantum Hall state with interaction driven 
edge canted antiferromagnetism. The local edge antiferro¬ 
magnetism, gapes out the original counter propagating spin 
filtered edge state. Band structure of a zigzag ribbon with 
spin orbit, whose edge magnetic moments driven by interac¬ 
tion are of plane (c) and in-plane (d). The in-plane magneti¬ 
zation destroys the gapless states due to breaking of TR and 
spin mixing, whereas the off-plane does not because Sz is still 
conserved. 


ing an AF ground state at zero in-plane fie ld, turning 
it into ferromagnetic above the critical field .1^^ In this 
situation, the edge does not play a key role, and the 
conductance only depends on the bulk magnetic order. 
A second scenario appears when cons idering that, ac¬ 
cording to theoretical predictions,the order in the 
bulk is actually ferromagnetic due to the long range 
Coulomb interaction. Starting in this situation, two 
counter-propagatin g sp in polarized states appear on the 
edge of the system.!^ Those states, can be highly vul¬ 
nerable to the local electron-electron interactions, and a 
local canted antiferromagnetic gap can be opened, gap¬ 
ing out the states and destroying the quantum Hall state. 
Thus, the critical field needed to reach the QSH is the 
one needed to recover ferromagnetism on the edge. In 
the context of the Hubbard model, it is found that at 
medium interactions (U = 2t), even in the presence of 
an exchange field, the bulk of the system can be in a fer¬ 
romagnetic state, at the same time the edge remains a 
canted antiferromagnet, as shown in Fig[^. 


VIII. EXPERIMENTAL EVIDENCE 

All the previous discussion was focused on our theo¬ 
retical understanding of a variety of different edge states 
in graphene and graphene like systems. Here we briefly 
touch upon the experimental evidence that has been ob¬ 
tained so far. 


A. Zigzag Edge states and edge magnetism 

We first discuss the observation of zigzag edge states 
and edge magnetism. Scanning tunneling microscopy 
(STM) had permitted early on to identity zigzag and 
armchair edges in graphitd^^. A fairly convincing 
evidenc^^ of the existence of localized edge states at chi¬ 
ral edges of ultrahigh quality graphene ribbons obtained 
by unzipping carbon nanotube^^was obtained by means 
of scanning tunneling microscopy. The same group also 
showed the ca pabi lity to engineer the edges by means of 
plasma etchinglSllof graphene nanoribbons (GNR). They 
found that hydrogen-plasma-etched GNRs are generally 
flat, free of structural reconstructions, and terminated by 
hydrogen atoms with no rehybridization of the outermost 
carbon edge atoms. Both zigzag and chiral edges show 
the presence of edge states. Edge states at the interface 
bet ween graphene and h-BN have also been imaged with 
STIVP^. 

The direct evidence for edge magnetism is less abun¬ 
dant, and most of the claims are based on indirect evi¬ 
dence. Part of the problem is that conventional magnetic 
probes, such as SQUID and EPR, have a resolution limit 
in the range of 10^^ fiB tha t wo uld require graphene flakes 
beyond current capabilitie j^^^l Therefore, such measure¬ 
ments are performed on solutions of graphene processed 
by various techniques, including unzipping carbon nan¬ 
otubes whose grow very often requires ferromagnetic cat¬ 
alyst. Therefore, the observation of r oom temperature 
ferromagnetism in graphene ribbonJ^^should be revised 
under the light of the recently worl^^SH discussing arti¬ 
facts for analogous claims for highly oriented pyrolitic 
graphiteP^. 

These shortcomings highlight the need to study edge 
magnetism in graphene using som e sort of local probe, 
preferably spin polarized STAf^i^J, or some type of nano- 
magnetometry technique, such as NV center detec tors, 
that can be used to probe local moments on a surface^^. 
A non-magnetic local probe was used by Tau and 
coworkerpl who compared their fairly extensive STM 
spectroscopy, both along the edge of the ribbon and in 
the direction perpendicular to the edge with mean field 
Hubbard model calculations where edge magnetism was 
included. The good agreement between theory and ex¬ 
periment is certainly a hint that edge magnetism was 
present. Further computational work using density func¬ 
tional calculations have confirmecP^ the presence of edge 
magnetism in hydrogenated graphene ribbons deposited 
on Au(lll). However, they also showed edge magnetism 
would not survive on other substrates, such as Gu(lll) 
and Ag(lll). 

The recent claims of room temperature edge ferromag¬ 
netism in graphene ribbond^ are also based on STM 
spectroscopy and comparison of the data to well estab¬ 
lished theoretical results. Unlike the work of Tau et aP^, 
a detailed evolution of the I{V) curves along the ribbons 
is not provided. In addition, edge magnetism is definitely 
not expected to survive up to room temperatur^^, al- 
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though the effect of long-range order depletion due to 
quantum and thermal fluctuations on the spectral func¬ 
tion of the quasiparticles has not been studied. A local 
probe of magnetism would be desirable in order to back 
up these extraordin ary claims. 

The recent reportP^^ of ballistic transport with a con- 

2 . . . 

ductance of G = ^ in epitaxial graphene nanoribbons 
would imply time reversal symmetry breaking compati¬ 
ble with edge magnetism. Great care was taken by the 
authors of this paper to rule out a great number of possi¬ 
ble experimental artifacts and the understanding of their 
data remains a very interesting challenge for theory. 

The finite spin of small triangulate molecules, with 
S' = i and S = 1 has been observed experimentally, in 
macroscopic samples, using electron paramagnetic reso¬ 
nance techniques. 


B. Edge transport in topological phases in 
graphene like systems 

The observation of the quantum Hall effect in graphene 
with perfect quantization of the Hall conductance pro¬ 
vides an extremely convincing, albeit indirect, evidence 
of the existence of Quantum Hall edge states. The Quan¬ 
tum Spin Hall predicted by Abanin et has been re¬ 
cently found experimentallj^^ by A. Young and cowork¬ 
ers. They verified, by means of local capacitance mea¬ 
surements, the opening of a band-gap in bulk and the 
edge nature of the transport by means of a floating 
gate techniqud^. The edge conductance in the QSH-like 
phase was G = I.SGq, below the ideal case of G = 2Go 
which suggest that some sort of spin-flip backscattering 
mechanism is present in the sample. 

IX. CONCLUSIONS AND OUTLOOK 

We have provided an overview of the electronic prop¬ 
erties of edge states in graphene, and related materials. 
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using as a guideline the standard tight-binding model for 
electrons in a honeycomb lattice. Remarkably, small ad¬ 
ditions on the same model that describes the main prop¬ 
erties of graphene, including the observed quantum Hall 
phase, have led to the prediction of the Quantum Spin 
Hall and Quantum Anomalous Hall phases, inspiring the 
search of graphene based topological phases. Therefore, 
small modifications of graphene, via proximity or func¬ 
tionalization, are expected to drive graphene into these 
exciting topological phases. 


Edge states play an important role in the study of some 
of the most fascinating electronic properties of graphene 
and graphene like materials. The study of graphene 
edge states has ramifications in many branches of modern 
material science, including graphene spintronics, organic 
chemistry, quantum topological phases, and metrology, 
just to mention a few. 

Combined with the discovery of new graphene-like ma¬ 
terials, and the construction of artificial structures where 
different two dimensional crystals are stacked and put in 
close proximity to materials with various types of elec¬ 
tronic order, such as ferromagnetism and superconduc¬ 
tivity, we can expect that the study of the electronic 
properties of edge states will bring many exciting results 
in the years to come. 
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